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ABSTRACT 

^ , Context. The small-scale distribution of dark matter in galactic halos is poorly known. Several 
•I— I . 

, studies suggest that it could be very clumpy, which turns out to be of paramount importance 
■ when investigating the annihilation signal from exotic particles (e.g. supersymmetric or 
Kaluza-Klein). 

Aims. In this paper we focus on the annihilation signal in positrons. We estimate the associated 
uncertainty, due to the fact that we do not know exactly how the clumps are distributed in the 
Galactic halo. 

Methods. To this aim, we perform a statistical study based on analytical computations, as well 
as numerical simulations. In particular, we study the average and variance of the annihilation 
signal over many Galactic halos having the same statistical properties. 

Results. We find that the so-called boost factor used by many authors should be handled with 
care, as i) it depends on energy and ii) it may be different for positrons, antiprotons and gamma 
rays, a fact which has not received any attention before. As an illustration, we use our results 
to discuss the positron spectrum measurements by the HEAT experiment. 
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1. Introduction 



Most observations of cosmological interest ean be accounted for by assuming that our 
Universe contains a large amount of non-baryonic matter, usually referred to as dark 
matter. Indeed, the mean density of matter f2,„ can be consistently estimated to be r2,„ ~ 
0.23 from many observations, whereas the baryonic density fif, inferred fr om primordial 



nucle osynthesis, from Cosmic Microwave Background (CMB) anisotropics ijSpergel et al 



200fi|) , Large Scale Structures and by direct observations of luminous matter is an order 
of magnitude lower, namely il,i,h^ « 0.0223. 

Among the several possible solutions to the dark matter problems, the hypothesis that 
it could be made of a weakly interacting fundamental particle of a new kind (hereafter 
wimp for weakly interactive massive particle) has received considerable attention. This 
is partly due to the fact that this hypothesis can be tested experimentally. In particular, 
the detection of the annihilation products of such exotic particles would be a great 
achievement, and an important fraction of the astroparticlc physics community is involved 
in that quest. 

However, assuming that wimps actually do exist (see e.g. 



Bertone et al 



12004) for 



a nice review on dark matter), their nature is basically unknown. Some constraints can 
be inferred from high precision cosmological observations such as the CMB, but several 
particle physics models provide candidates whose properties are consistent with these 
observations. Extensions of the standard model of particle physics, such as supersymme- 
try and Kaluza-Klein theories, naturally offer such candidates. The lack of information 
about the nature of the wimp may translate into a very large uncertainty on its spatial 
distribution. There are two main physical reasons for that. 

First, the gravitational collapse of primordial density inhomogeneities that leads to 
the formation of cosmic structures is characterized by a small scale cut-off, due to several 
physical effects. To begin with, the particles enduring the collapse may interact elasti- 
cally with other species or between themselves, which is responsible for diffusion. Then, 
after their interactions become negligible, they become free to move out of the collaps- 
ing region: this i s free- s treaming. A general discussion of these effects can be found in 

may st r ongly depend on the n ature and 



Berezinskv et al 



(2003). The res ulting cut-off 



properties of the wimp (se e e.g 



Boehm et al 



stance, the recent study by 
3 X 10-3 Mm to 3 X 10-1 M, 



Profiimo et al.l 



l200f|) 



(200l|); 



Hofmann et al 



200l|)). For in- 



givcs protohalo masses ranging from 



0- 
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Then, the structures evolve, merge and can be partially disrupted by tidal forces, so 
that the current cut-ofF in the spectrum of clump masses corresponds to t he smallest sur- 
viving clumps. The situation is still unclear, as numerical simulations bv lDiemand et al 



(20051) show ed that clumps as small as 10 ^ M0 could survive disruption, while analyt- 



ical work by 



Berezinskv et al 



showed that structures smaller than 10 Mrr, were 



disrupted. The possibility th at tidal interaction with stars may play an important r ole 



has also been hotly debated ijZhao et al 



2005a: 



Moore et al 



2005 



Zhao et al 



2005b) 



Different very competent experts of the domain provide very different descriptions of 
the dumpiness of galactic halos, and the aim of this paper is not to tackle with this issue. 
Instead, we consider the wide range of possibilities as the starting point of our analysis. 
The quantitative amount of dumpiness is of paramount importance as it enhances the 
annihilation rate of wimps and increases the detection prospects. 

In most studies, dumpiness is taken into account by a general, energy-independent 
multiplicative number called boost factor, by which the signal computed from a smooth 
dark matter distribution should be multiplied. This is not correct and we show in this 
paper that the effects of dumpiness cannot be described by such a unique number. 
Moreover, this is a stochastic problem, in the following sense: general hypotheses about 
the statistical properties of the distribution of clumps in the Galactic halo can be made, 
but the exact position of every clump is unknown. In some cases, the expected signal 
from a given type of wimp can be quite sensitive to the precise position of the Earth 
relative to the nearest clumps. 

The aim of this paper is to study the effect of the halo dumpiness on the annihilation 
signal, focusing on the case of positrons. Taking advantage of analytical computations 
and numerical simulations, we investigate the statistical properties of the annihilation 
signal. We show that, at variance with the assumptions of most studies, the dumpiness 
factor depends on energy and is not the same for positrons as for gamma rays. We also 
show that even if the average properties (averaging being meant over a large number of 
realizations of our Galactic halo) of a clumpy halo may be well described by the usual 
boost factor, the deviations from this average may be very large and the ability to predict 
a signal from a model may be consequently deteriorated. 

The importance of dumpiness in deter mining the dark matt er annihilation signal 



in positrons has alread y be en assessed by 



Hooper et al 



Baltz fc Edsiol l|1999|) . and further studied 



( 2004^ and iHooper & Kribsl (|2004 . The possibility that the positron 
excess observed by HEAT could be due to a single nearby clump had been raised. The 
probability of such a sit uation was estimated to b e low (about 10^^). More recently. 



this proposal resurfaces (jGumberbatch fc Silk 



20061) . As an illustration of the methods 
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developed in this work, wc show that the odds for such an occurrence are even lower 
than Hooper et al's estimations. 



2. The effective boost factor 

For the sake of dcfiniteness, we first consider the case of wimps annihilating into positrons 
and electrons at a given energy - the source sp ectrum of positrons c an be considered mo- 
nocncrgctic. In Kaluza-Klcin inspired models ( Servant fc Taitl2003 L dark matter species 
may substantially annihilate into electron-positron pairs with a branching ratio as large 
as ~ 20%. The positron production rate Pe+ counts the number of annihilations taking 
place per unit volume at some point x 

P^+{x)=6{a,n,Jxx-e+e-)v) {P^} (1) 



where the S term is equal to 1/2 for Majorana particle, taking into account the fact that 
these particles are not discernible, whereas it is equal to 1/4 for Dirac particles, taking 
into account the fact that the density of particles and antiparticles is p/2 and not p. The 
contribution of the infinitesimal volume d^x located at point x to the flux at the Earth 
- in units of cm^^ s^^ sr^^ GeV~^ - of the resulting positrons with energy E may be 
expressed as 

dcf) = SGe+ {xq,E ^ x,Es) 1^^^-^^ d^x , (2) 

where the quantity S depends on the mass density of reference po and on the specific 
features of the high energy physics model at stake 

5=^«e+(i^)(^an„(xX-e+e-) t;) (^1 . (3) 
47r ^ '1^ J 

The velocity of the positron with energy E at the Earth is denoted by We+ . The Green 
function Ge+ is discussed in section |3 It describes the probability that a positron pro- 
duced at point X with energy Eg reaches the Earth with a degraded energy E. As wimps 
are at rest with respect to the Milky Way, the energy Eg is equal to the parent parti- 
cle mass m^. At this stage, we would like to keep our discussion as general as possible. 
Because our formalism should easily be extended to any charged species - to antiprotons 
or antideuterons for instance - the positron propagator will be denoted more simply as 
G {x, E). The total positron flux at the Earth results from the integral over the galactic 
DM mass distribution p{x) 

(f>^S f G{x,E) 1^1' d'^x . (4) 

Should the DM halo be smoothly distributed with mass density ps, the positron flux 
would be given by relation where the wimp distribution is now described by ps 

<Ps=S I G{x,E) [P^V d^x . (5) 

JDMhalo I Po J 
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In the literature, the effects of dumpiness have been so far accounted for by shifting 
upwards the flux cjjg. The multiplicative factor is called the boost. It acts as a constant 
of renormalization by which the flux 0^ generated by a smooth DM halo should be 
multiplied in order to take into account the enhancement of the wimp annihilation rate 
inside substructures. That procedure has been widely used in the past but is shown to 
be wrong in the present paper. In the following we discuss the method that must be 
followed in order to compute correctly the signal (p at the Earth. 

We assume that substructures ~ whose density profile inside the i-th clump is 5pi{x) 
- float inside a smoother background with mass density p' ^ which is a priori different 
from ps introduced above. The halo density p can be written as 



P's + ^5p, , (6) 



and each clump has a mass 



M, = / d^'x Sp^ (x) . (7) 

J ith clump 

Because wimp annihilation involves the square of the DM mass density, the production 
of positrons inside the i-th protohalo is enhanced with respect to the situation where 
that substructure would be diluted in the surroimding medium. Should the latter be 
homogeneously spread with a mass density ph (which will correspond to fps below), the 
enhancement would be given by the boost factor Bi which we define as 

d^x 5p} {x) =M,x Ph B, . (8) 

ith clump 

That relation does not mean that the annihilation signal scales linearly with the clump 
mass. The boost factor Bi actually takes into account the inner DM distribution so that 
various profiles for 5 pi can lead to very different values for Bi. The rel evant quantity 



turns out to be the effective volume Bi Mi/ ph- In the case of model (B) of (jBertone et al 



2005|) where the DM clumps have been accreted around intermediate-mass black holes, 
the average value for that crucial factor is ~ 4 x 10^ kpc'^ even if the spike radius is only 
~ 1 pc. Relation © translates into the positron flux at the Earth 

= + , (9) 

whose component 



is produced by the constellation of DM protohalos that pervade the Milky Way. The 
signal (f) at the Earth is therefore enhanced by a factor of B = <f>/(l)s with respect to the 
situation where the DM halo is completely smooth with mass density Ps- Many clump 
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distributions arc possible and lead to different values for the boost B. The distribution 
inside which we are embedded is of course unique. Unfortunately, we know little about it. 
In order to predict the set of plausible values for the boost B, we are forced to consider the 
vast ensemble of all the possible DM substructure configurations. Our lack of knowledge 
limits us to merely derive trends for the boost. The analysis of how B is statistically 
distributed is postponed to section [4.11 Instead we now focus on its average value Bas 
which suffices when its variance is small. To proceed further, a few simplifications are 
nevertheless helpful. 

(i) We will first assume that clumps are practically point-like. This hypothesis is 
expected to be valid when the propagation distance is large compared to the size of the 
clump. As the volume of the galaxy filled by the DM substructures becomes negligible, 
the halo density p becomes 



P . 



(11) 



and the smooth component 0'^ of the flux is given by relation © where the mass density 
Ps is now replaced by p'^. Moreover, the positron flux ipi which the clump located at 
position Xi yields, simplifles into 



PO 



(12) 



where d = G{xi,E). 

(ii) The boost factor Bi at the source should vary from one protohalo to another even 
if the mass Mi is assumed constant. To commence, the inner regions of the Milky Way 
have presumably collapsed earlier than its outskirts, dragging with them substructures 
whose concentrations are stronger than for the galactic periphery. We could expect to 
have larger values of Bi inside the solar circle. However, clumps that mo ve near the galac- 



tic ce nter experience strong tides that could significantly reshape them ([Berezinskv et al 



2003|) . It is not unconceivable that clumps partially evaporate like globular clusters which 
exhibit characteristic tidal tails. If that effect is dominant, the clump mass is reduced 
and probably the boost factor too - if the density profile of the substructure readjusts 
itself accordingly. It is therefore quite difficult to predict how the clump boost factor Bi 
varies with position. In order to simplify the discussion, we assume that all the clumps 
have the same mass Mi = Mc and the same boost factor Bi = B^. The fi rst hypothesis 
is supported by numerical simulations that indicate I Diemand et alJLoOfi ) that the mass 
function of substructures is a self-similar power law of slope dn{M)/d\ogM oc M^^ 
and is actuall y dominated by th e lightest clumps. The latter hypothes i s is a priori more 



questionable l|Zhao et al 



2005a: 



Moore et al 



2005 



Zhao et al 



2005b : 



Berezinskv et al 
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2003). It is nevertheless a reasonable choice insofar as the effects of a granular DM dis- 
tribution on the flux of positrons will be shown to be mostly local. The actual value of 
Bi should not vary much in the solar neighbourhood and we can safely consider that it 
is constant. The expression for the flux (f) simplifies into 

BcMc 



Po 



Eg. 



(13) 



(iii) A fraction / of the total DM halo is in the form of substructures embedded 
inside a s mooth component wit h mass density p'^,. In the intermediate-mass black hole 



scenario of 



Bertone et al 



in the IjDiemand et al 



(12005^ . the fraction / is so small that p' , ~ ps- On the contrary, 
boOli ) simulations, a value as large as / ~ 0.5 is found with 
a preponderan ce of small-scale c lumps which should trace the smooth DM density as 



ascertained in 



Berezinskv et al 



1 2003f) . The mass density p' ^ could be quite different 
from Ps but its contribution (p'^ to the overall signal (j) is small. We will assume for 
simplicity that 

p'., = {l-f)Ps , (14) 

where / is constant all over the Milky Way. The corresponding flux ratio (t)'s/4>s ^ which 
should not exceed 1 in any case - is now given by the factor (1 — /)^. 

(iv) A number Nh of DM substructures pervade the Milky Way halo. In this analysis, 
we will not consider the fluctuations of that number. The probability that one of those 
lies at point x is controlled by the distribution p{x). The number of clumps which the 
volume d?x contains on average is 



{dn) = Nhp{x) d^x . 
We infer an average flux at the Earth 

(0> = (l-/)^0. + 5^ l^G^ , 

where the average sum over the Green functions is given by the integral 



Eg. 



a {x, E) (dn) 



(15) 



(16) 



(17) 



'DM halo 

For illustration purposes, we have chosen in our numerical examples a particular clump 



distribution. Inspired by (|Diemand et al 



2005|) . we have assumed that protohalos trace 



the smooth distribution of dark matter with 



p{x) 



Ps (x) 
Mh 



(18) 



where Mh is the mass of the DM Milky Way halo. We stress that our analysis does 
not depend on that spcciflc choice and is completely general. Considering a different 
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distribution p (x) - with no relation to the mass density in particular - would not 
qualitatively affect the main conclusions of our analysis. 

Keeping this in mind, we can proceed with our illustrative choice for p [x) and derive 
the effective boost 

Beff(i?) ^ ^ = (1 - + /Bc^ , (19) 

where the integral 2„ is defined by 



In{E) = f Gix,E) 

JDMhalo 



Clump Fraction f = 0,2 
Isothermal DM Halo with a core 




Positron Energy E at the Earth [GeV] 
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10 
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1 

[GeV] 



(20) 



Fig. 1. The effective boost factor iJoff is featured as a function of the positron energy E 
in the case of a 100 GeV line. A fraction / = 0.2 of the DM distribution is in the form 
of substructures whose individual boost factor - relative to the solar neighbourhood 
density - has been varied from 3 to 100. An isothermal halo - panel a - and a NFW 
profile - panel b - are considered. They illustrate the influence of the central profile 
index. The increase of i?off is noticeable especially around i? 10 GeV. 

Although the boost at the source is fixed, the boost of the signal at the Earth 
Bcff depends on both the nature and the energy of the cosmic ray species through the 
Green function G and the integrals X\ and Xi . We reach the conclusion that as the flux 
0s is shifted upwards as a result of dumpiness, it also experiences a spectral distorsion 
insofar as i?off is energy dependent. This property has never been noticed before and is 
illustrated in the two panels of Fig.nwhere the case of a 100 GeV positron line is featured. 
A fraction / = 0.2 of the DM halo h as collapsed in the for m of clumps whose boost factor 
Be is varied from 3 to 100. In the jpiemand et allboOfJ) numerical simulations, such a 
value for the fraction would correspond to a minimum mass scale of 10"* Mq. The latter 
lies typically at the lower tip of the range of protohalo masses which we have used in 
our examples. As for the boost S^, the values quoted in the literature vary from a few 
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1 Berezinskv et al 



2005). The 



2006|) up to over two orders of magnitude (jDiemand et al. 
mass density of reference po ti^^s been set equal to the solar neighbourhood value of 
Ps{Q) — 0.3 GcV cm~'^. The increase of the effective boost factor with positron energy 
is clear in both panels. Near the line - in the region where E tends to the input energy 
Es - the positron Green function G probes only a small region of the Milky Way halo 
around the solar system. With our definition of pQ, the integral ratio I1/I2 boils down 
to unity and Bcs may be approximated by ~ / -Be- If £■ is now varied from its upper 
limit £^5 downwards, larger portions of the halo come into play in the integrals Xi and X2 
with the result of decreasing their ratio. That effect is quite obvious in panel a) where an 
isothermal profile is assumed with core radius Ocoro = 0.5 kpc. The DM density ps within 
1 kpc from the galactic center is larger than in the case of a NFW distribution and the 
relative increase of X2 " where the square of ps is relevant and not merely ps alone - with 
respect to Xi is more pronounced. It is possible that the energy dependence of Bcft which 
we have discovered could strengthen the case of DM particles as a plausible explanation 
of the still putative positron excess reported by HEAT. In the example featured in Fig.^ 
the largest spectral distortion is actually obtained for a positron energy E 10 GeV. 

Note that this distorsion effect should not be present in the case of gamma rays, 
whose propagation does not depend on energy. As regards antiprotons, the Green function 
already probes a significant portion of the DM halo, and wc anticipate a mild dependence 
of the boost factor on the energy. 

3. The positron propagator 

The departures of the positron flux from (f>s are expected to be large when the positron 
energy E is close to the production value Es- In this regime, the particles cannot have 
been produced far away. They mostly originate from a region close to the solar system 
inside which the distribution of clumps may significantly fluctuate. That is why wc have 
focused our analysis on cosmic ray positrons whose propagation throughout the galaxy 
is now briefly sketched. 

The master equation for positron propagation is the continuity relation 

d^J^ + OeJ^ = Q , (21) 

where Q denotes the production rate of positrons per unit of volume and energy. The 
space-time vector current is defined as 

- % ■ <^^' 

The time-component J" = dn/dE = ip {x, E) denotes the number density of particles per 
unit of volume and energy. The space current accounts for the scattering of cosmic rays 



10 Lavalle et al: Positron annihilation in a clumpy halo 

upon the inhomogeneities of the galactic magnetic fields which is described as a diffusion 
process with 

J = -K{x,E) W . (23) 
The energy component depends on the energy loss rate b through 



[^E)^-biE)} . (24) 



Above a few GeV, positron energy losses are dominated by synchrotron radiation in the 
galactic magnetic fields and by inverse Compton scattering on stellar light and on CMB 
photons. The energy loss rate b depends on the positron energy E through 

biE) = -f- . (25) 

We have set the energy of reference Eq to 1 GeV and the typical energy loss time is 
te = 10^^ s. The master equation (|21|) may be expanded into 

^ -V-{K {x, E) Vi,} - ^ {b{E) i,}^Q {x, E) . (26) 

In order to simplify the discussion, steady state is assumed and the space diffusion coef- 
ficient K is taken to be homogeneous with the energy dependence 

X(i?) = /lo 1^1" . (27) 

The diffusion coefficient at 1 GeV is Kq = 3 x 10^^ cm^ s~^ with a spectral index of 
a = 0.6. The master equation (|26|) simplifies into 

Aoe"A7^ + |-( — + Q = , (28) 
oe {te J 

where e denotes the ratio E/Eq- 



Equation H28|) may be solved with the 



Baltz fc Edsiol l|1999|) clever trick which consists 



in translating the energy E into the pseudo-time 



tiE)=TE i^viE) = ^1 . (29) 

The energy losses which positrons experience boil down to a mere evolution in this 
pseudo-time so that the propagation equation H28|) greatly simplifies into 

^ - K„Ai, = Q{x,i) . (30) 
ot 

The space and energy positron density is now ip ^ ip whereas the positron production 
rate has become Q = gS-a q fsfQ^j^-g i^j^^t both tf; and Q have the same dimensions as 
before because e is dimensionless. Without any space boundary condition, equation (|30|) 
may be readily solved. If a drop is deposited at the origin of the coordinates at pseudo- 
time is = 

Q{xs,is)=S''{xs)S{is) , (31) 
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the subsequent diffusion in an infinite 3D space would result into the density at position 
X and pseudo-time t given by the well-known Green function 

=G(a;,i^ 0,0) =0(t) {4^X0?}"'^' cxp |-^^| , (32) 

where r = \x\. The general solution of equation H3U|I may be expressed with the Green 
function G as 

i,{x,t)^ f' dts [ (fxsG{xJ^xs,ts) Q{xs,is) , (33) 

Jts=0 J 

and translates into 

^{x,E)= dEs d^xsG,+ {x,E^xs,Es)Q{xs,Es) . (34) 

JEs=E J 

The positron propagator may be obtained from G through 

G,+ {x,E^xs,Es) = ^G{x,i^xs,is) , (35) 

where the connection between the energy E and the pseudo-time t is given by rela- 
tion (|29|l . In the case of monochromatic positrons, the production rate is 

Qix,E) = P,^{x)d{E-Es) , (36) 

and the positron space and energy density at the Earth may be expressed as 

^ixQ,E) ^eiEs- E) J d''xsGe+ixQ,E<-xs,Es) Pe+ixs) ■ (37) 

Equation (0) is based on this relation. 

The diffusive halo inside which cosmic rays propagate before escaping into the inter- 
galactic medium is pictured as a flat cylinder with radius i?gai ~ 20 kpc and extends along 
the vertical direction from z = —L up to 2: = +L. We have assumed here a half-thickness 
of L = 3 kpc. Without any boundary condition, the propagator G would be given by 
the 3D relation (|32|) . However, cosmic rays may escape outside the diffusive halo and G 
should account for that leakage. In spite of the boundary at i?gai, we have assumed that 
cosmic ray diffusion is not limited along the radial direction but takes place inside an infi- 
nite horizontal slab with thickness 2L. We have nevertheless disregarded sources located 
at a radial distance R larger than i?gai. Indeed, because their energy is rapidly degraded 
as they propagate, positrons are observed close to where they are pro duced. Our radial 



treat ment is justified because positrons do not originate from far away l|Maurin fc Taillet 



2003^ . Even in the case of antiprotons for which the galactic propagation range is signifi- 



cantly larger than for positrons, the effects of the radial boundary down at the Earth are 
not significant insofar as cosmic ray species tend to leak above and beneath the diffusive 
halo at z = ±L instead of traveling a long distance along the galactic disk. The infinite 



12 Lavalle et al: Positron annihilation in a clumpy halo 

slab hypothesis allows the radial and vertical directions to be disentangled in the reduced 
propagator G which may now be expressed as 

G {x, i ^ xs, is) = , ^^l^ - exp (- 1 V {z, i ^ zs, is) , (38) 

where f = t—ts- The radial distance between the source xs and the point x of observation 
is defined as 

1 /2 

R= {{x-xs)^ + iy^ys)^} ■ (39) 

Should propagation be free along the vertical direction, the propagator V would be given 
by the ID solution Vm to the diffusion equation (|3l)|l 

V (z, i ^ zs, is) = Vizj {z, i <- zs,is) = -7== exp | - ~ 1 . (40) 

But the vertical boundary conditions definitely need to be implemented. Wherever the 
source inside the slab, the positron density vanishes at z = ±L. A first ap proach relies on 
the m ethod of the so-called electrical images and has been implemented by 



Baltz fc Edsic 



(1999). Any point- like source inside the slab is associated to the infinite series of its 
multiple images through the boundaries at z = ±L which act as mirrors. The n-th image 
is located at 

z„ = 2L7i + (-1)" zs , (41) 

and has a positive or negative contribution depending on whether n is an even or odd 
number. When the diffusion time f is small, the ID solution (|40|l is a quite good approx- 
imation. The relevant parameter is actually 

C - , (42) 

4 Aqt 

and in the regime where it is much larger than 1, the propagation is insensitive to the 
vertical boundaries. On the contrary, when C is much smaller than 1, a large number of 
images need to be taken into account in the sum 

V{z,i^ zs,is) ^ (-1)" ViD (z,?^ z„,ts) , (43) 

n — — OO 

and convergence may be a problem. It is fortunate that a quite different approach is 
possible in that case. The ID diffusion equation H30() actually looks like the Schrodinger 
equation - in imaginary time - that accounts for the behaviour of a particle inside an 
infinitely deep ID potential well that extends from z = —L to z — +L. The eigcnfunctions 
of the associated Hamiltonian are both even 

(^„(z) sin{fc„ (L- |z|)} (44) 

and odd 

(^;(z) = sin{/c;(L-z)} (45) 
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functions of the vertical coordinate z. The wave- vectors /c„ and /c^ are respectively defined 
as 

/ 1 \ TT 7r 

\ Z J Lj Lj 

The vertical propagator may be expanded as the series 
V{zJ^zsJs) ^Y. J {cr^''^^nizs)ipni^) + e-\^^ip'^izs)ipU^)} , (47) 

n — 1 

where the time constants A„ and are respectively equal to -ftTg a-nd k'^^ . In the 
regime where C is much smaller than 1 - for very large values of the diffusion time f - 
just a few eigenfunctions need to be considered for the sum l|47|l to converge. 



4. An analytic approach of the cosmic ray flux fluctuations 

4.1. The random flux (pr and its variance 

The cosmic ray flux at the Earth contains the random component 

V^-S^Y.^^ , (48) 

which is produced by the constellation of DM clumps inside the Milky Way halo. Before 
embarking into our discussion, a few remarks are in order. 

(i) The actual distribution of DM substructures is of course unique and so is the 
cosmic ray flux which it generates at the Earth. We will nevertheless consider it as one 
particular realization among an essentially infinite number of different other possible 
realizations. We furthermore assume that clumps are distributed at random and that the 
set of all their possible distributions makes up the statistical ensemble which we consider 
in this section. The aim of our analysis is to investigate how strongly the flux 0^ may 
fluctuate as a result of the random nature of the wimp clump distribution. We will derive 
the associated cosmic-ray flux variance ar defined as 

= (02) _ (0^)2 . (49) 

The variance <7r turns out to be an essential tool. Because the total flux (/) and its random 
component 0,. are shifted with respect to each other by the constant quantity 

^-<f>r = (l-ff , (50) 

the variance of the former is given by 

al = (02) - (0)2 = (02) _ (^^)2 ^ ^2 ^ (51) 
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The effective boost itself B^s which has been discussed in section |21 is an average value 
around which the true flux enhancement B = 6/6^ fluctuates with the variance 



'^B^ ^ = :r ■ (52) 

<Ps <Ps 

Therefore, the determination of Ur leads immediately to the boost fluctuations as- 

(ii) We will furthermore assume that clumps are distributed independently of each 
other. The problem is then greatly simplified because we just need to determine how 
a single clump is distributed inside the galactic halo in order to derive the statistical 
properties of an entire constellation of such substructures. In particular, the average 
value {(j)r) of the random component of the cosmic ray flux is readily obtained from the 
average flux {ip) produced by a single clump through the relation 

, (53) 

where Nu denotes the total number of clumps to be considered. The variance cr^ - which 
is the crucial quantity as regards the flux fluctuations - may also be expressed as 

al^Nua^ ^Nh{{v'') - {V?) ■ (54) 

(iii) The set of the random distributions of one single clump inside the domain 'Dh 
forms the statistical ensemble T which we need to consider. An event from that en- 
semble consists in a clump located at position x within the elementary volume (fix. Its 
probability dP will be assumed to follow the smoothed DM mass distribution ps so that 

dP = p{x) d^x^^^d^x . (55) 

The domain T>h over which our statistical analysis is performed is so large that the total 
number Nh of clumps which it contains is essentially infinite. That region Vh behaves 
therefore like a so-called thermostat in statistical mechanics. It encompasses of course 
the diffusive halo and may even be much bigger. It may be thought - but not exclusively 
- as the entire Milky Way DM halo. Its actual size has no importance because it will 
disappear from the final results in the limit where the ratio 1/Nh is negligible. The only 
requirement is that Nh should be much larger than the typical number Ns of clumps 
that effectively contribute to the signal (pr at the Earth. The domain Vh contains the 
total DM mass Mh ~ & fraction / of which consists in Nh identical clumps so that 

NHMc = fMH . (56) 



We are now ready to derive the probability distribution V{ip) associated to the signal 
if which a single clump generates. The statistical properties of the random variable (p {T} 
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translate those of the statistical ensemble T itself. More precisely, the probability function 
Viip) is related to the space distribution p (x) through 

V{ip)dip^dP= [ p{x) (fix . (57) 

The subdomain TD^p over which the space distribution p [x) should be integrated in the 
previous expression yields a flux at the Earth comprised between Lp and (f -\- dip {Dh is 
thus the union of all 'Dip). In the case of positrons, the probability distribution V{'p}) will 
be shown in section [4.21 to concentrate around a flux ip equal to 0. The average value - 
over the statistical ensemble T - of any function J- that depends on the flux ip may be 
expressed as 

{T) = I T{p) V{p)dp= [ T {p {x)} p {x) d^x . (58) 

J-Dh 



In particular, the flux which a single clump yields on average at the Earth is readily 
derived from the integral 



{p) = p{x) p{x) d 



where X„ has been defined in relation H2U|I . The average value of the random flux (j)^ 
implies Nh clumps and expression l|53|l - with the help of relation H5f)|l - leads to 

((br) S J Bell II , , 

^ = . =!Bc^ , (60) 

and to formula (|19|) . 

Starting from the definition H54|). the variance tr^ may be derived in the same spirit 
with 



{4>r? Nh i^)^ Nh 
With the help of relation H58() , the mean square of the single clump flux may be expressed 
as 

{^')= [ p'ix)p{x)d^x=^^^&Ji, (62) 
Jvh Po Mh 

where the integral is defined as 

Jn{E) = f {x,E) leiM.y d^x . (63) 

A straightforward algebra leads to the relative variance 

(0.)' P^Nh Tl Nh ^ fpo If ' ^ ' 

Because the domain Vh is so large - remember that both Vh and the Milky Way DM 
halo encompasses the diffusive halo - we can safely drop the ratio 1/Nh in the previous 
expression. 
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Positron Energy E at the Earth [GeV] 



Fig. 2. The relative variance (T,./((/),.) of the random component of the positron flux at 
the Earth - solid lines - and its hard-sphere approximation - long-dashed curves - are 
featured as a function of the positron energy E for three different values of the clump 
mass Mc- The injected positron energy Es has been set equal to 100 GeV. A NFW profile 
with typical scale 25 kpc has been assumed. At fixed clump mass, the variance increases 
with E and matches its hard-sphere approximation above ^ 40 GeV. As the number of 
clumps is decreased, the curves are shifted upwards by a factor of 1/\/Nh oc \fWc. The 
relative variance (Jb/BoS of the boost factor is also displayed ~ short-dashed curve. In 
the limit where the clump boost factor Be is large - a value of 100 has been assumed 
here - as /Beg and Gr/ {4>r) are approximately equal. 

The positron propagator of section |3| has been used in relation 1)64(1 in order to derive 
the solid curves of Fig. |21 At fixed Nh , the clump mass Mc is determined by equation 
1561 and the relative variance Ur/ {4>r) mcreases with the positron energy E at the Earth. 
This behaviour will be explained in section with the hard-sphere approximation. The 
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Fig. 3. The effective boost factor B^s ^ black line - is plotted as a function of the positron 
energy E for an injected energy Es = 100 GeV. The 1-cr range of its fluctuations extends 
from -Bmin = Soft ^ cs up to -Bmax = ^off + cr B ■ At fixcd clump mass, that range opens 
up as E approaches the injected energy Es = 100 GeV. It also widens significantly at 
fixed positron energy E when the number of clumps is decreased. 



ratio cfr/{4'r) is proportional to 1/ \/Nh oc \/Mc, and weighted by an effective volume 
J'l/Xf. The curves are therefore shifted upwards when the clump mass is increased. The 
relative variance as/Bcs of the flux enhancement B = (j)/(j)s is also presented in Fig. |21 
In the limit where the individual clump boost factor Be is large - we have selected here 
a value of Be = 100 - the random component 0,- of the positron flux dominates over its 
smooth counterpart (1 — /)^ <f)s so that 
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That is why the solid lines and short-dashed curves of Fig.|21are quite similar. In Fig.|31 
the black central curve features the effective boost factor B^tf of a NFW halo and corre- 
sponds to the case Be ~ 100 of the panel b of Fig. ^ from which it has been extracted. 
The l-cr range of its fluctuations extends from i?min = Scff^crs up to Bmax = ^off + 
At fixed clump mass, that range opens up as E approaches the injected energy Es ~ 100 
GcV. The fluctuations in the positron signal increase significantly just below the positron 
line. The boost variance ctb is also proportional to l/y^/Vc oc VMc- That is why the fluc- 
tuation band broadens up as the clump mass is increased from lO"' up to 10^ M0. 

4.2. The flux distribution V{ip) of a single clump 

The positron flux at energy E < Es which a single clump located at position x generates 
at the Earth implies the propagator discussed in section 13 

ip{x) = S^^Ge+{xQ,E^x,Es) . (66) 
Po 

and may be expressed with the reduced Green function G as 

(fix) =S ^"^^^ G {xQ,i^ x,is) ■ (67) 

Po Eoe^ ^ 

When the substructure is very close to the Earth, the flux (p reaches a maximal value 
</5max that depends both on the clump properties through the effective volume BcMc/ po 
and on the specific features assumed for the DM particle through the factor S. Without 
any loss of generality, we can significantly simplify the discussion by considering the ratio 

= (68) 

finax Umax 

instead of the fiux ip itself. Wc therefore would like to derive the density of probability 
7'(<I>) associated to the reduced flux <I> as it varies from to 1. 

In Fig. 21 that distribution is presented for three typical values of the positron energy 
E at the Earth. The energy Es of the positron line has been set equal to 100 GeV and 
a NFW DM halo has been assumed. The solid curves correspond to the fully numerical 
calculation of V{^) based on relation (|57|l . The domain Vh over which the probability 
is normalized to unity is the Milky Way DM halo up to a radius of 20 kpc. That domain 
encompasses the diffusive halo outside which the cosmic ray density vanishes. Most of the 
probability is therefore contained in the low flux region and the density V{^) basically 
diverges at $ = 0. As the energy E increases towards £'5, the region of the diffusive 
halo that is probed by the positron propagator shrinks. That region corresponds to large 
values of the positron flux <&. As its volume decreases when E approaches Es, fewer 
clumps are involved in the signal and the corresponding probability lessens. Notice in 
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0,2 0,4 0.6 0,8 



Reduced flux ^ = (p/<Pme.x ^^'^^ a single clump 



Fig. 4. The density of probability V{^) is plotted as a function of the reduced flux 
$ — (p/(pma.x which a single clump generates. A NFW halo has been assumed with a 
scale radius of 25 kpc. The domain Vh over which the probability is normalized to unity 
is the Milky Way DM halo up to a radius of 20 kpc. The injection energy is £'5 = 100 
GeV. The smaller the positron energy E at the Earth, the larger the probability density 
for a non-vanishing flux. The fully numerical calculations - solid curves - are compared 
to the infinite 3D approximation 1)71(1 that corresponds to the long-dashed lines. 

Fig. 0] how the probability density drops when E is increased from 1 GeV up to 

50 GeV. The lower curve is reproduced in Fig. Altogether with the distributions $7'(<f>) 
and <i>^ V{^) whose integrals from $ = up to $ = 1 are respectively related to {(p) and 

When the positron energy E is close to the energy Es, the pseudo-time difference 
T ^ t — ts is so small that the diffusion is no longer sensitive to the vertical boundaries at 
z = —L and z = +L. The Green function G can be safely approximated by the gaussian 
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Fig. 5. The density of probability 7'(<i>) as well as and <i>^7'(<i>) are featured as 

a function of the reduced flux $ — (p/(fiinax for a positron energy at the Earth of 50 GcV. 

function (sec OH) 

G{xQ,i.-x,is) ={i7rKof}-^^^ cxpf^^j^^ . (69) 

This regime corresponds to large values of the parameter C - defined in relation (|42() - or 
alternatively to small values of the positron diffusion length Ad = V 4 Kq f . The latter is 
featured in Fig. as a function of E for three different values of the energy at source. In 
the case where Es = 100 GcV, the diffusion length Ad exceeds the thickness L below an 
energy of ~ 8 GeV. Above that limit, positron propagation is not affected by the vertical 
boundaries and the infinite 3D approximation (|69|) applies with a reduced flux $ that 
only depends on the distance r of the clump to the Earth 



<I>-cxp (-rVA^) 



(70) 
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Fig. 6. The positron diffusion length Ad decreases as the energy E at the Earth ap- 
proaches the energy Es of the Hnc. The long-dashed horizontal line corresponds to a 
diffusion length Ad equal to the thickness L of the diffusion layers. Below that limit, 
positron propagation is not sensitive to the vertical boundaries and the infinite 3D ap- 
proximation is valid. This regime corresponds to large values of the parameter ^ - see 
the definition (g^J- 



An analytic density of probability may be derived in that case with 



-P($) = 2 7rA^ 



(71) 



That relation corresponds to the long-dashed curves of Fig. 0] where a value of Mh ~ 
1.357 X 10^^ Mq has been found for the mass contained in the inner 20 kpc of the 
Milky Way DM halo. When the positron diffusion length Ad is smaller than the slab 
thickness L, relation H71|l is an excellent approximation to the density of probability 
V{^). As an illustration, we find a value of Ad = 1-26 kpc well below L = 3 kpc when 
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the positron energy E is equal to 50 GeV. It is no surprise therefore if the solid and long- 
dashed red lines of Fig. 0| are so well superimposed on each other. As E decreases, the 
diffusion length Ad becomes more and more sizeable with respect to L and the infinite 
3D propagator H69|l tends to overestimate the region from which the signal originates 
as well as the corresponding density of probability Notice how the long-dashed 

approximation lines are shifted upwards with respect to the solid true numerical curves 
in Fig. 01 As E decreases, the approximation ()71|l worsens and the disagreement with the 
correct result becomes more pronounced. 



5. The hard-sphere approximation 

In the limit where the infinite 3D approximation applies - actually for a large range of 
values of the positron energy E at the Earth - we can simplify further the propagator 
Ge+ and substitute the step function 

G {xQ,i <- ^r^Js) ^ (72) 

for the gaussian form Ht)9|l . The distance between the clump and the Earth is denoted 
by r = \x — Xq\. According to this hard-sphere approximation, the Green function G 
reaches the constant value 1/Vs inside the sphere Vs of radius rs and volume Vs — 
whose center coincides with the Earth - and vanishes elsewhere. Both expressions (|69(l 
and (ITU are normalized to unity. The integral over the full 3D space of the square of 
those Green functions should also be the same. This condition translates into 

1-= j &<fx , (73) 

and leads to the volume 

Vs = (\/2^Ad)^ . (74) 

In spite of its crudeness, the hard-sphere approximation turns out to be quite powerful 
and is an excellent tool to understand the salient features of the statistical properties of 
the clump distribution and of its flux. The associated density of probability has little to 
do with the curves of Fig. 01 or with relation ^7\\ . It is actually a bimodal distribution 
with 

V{>^)=p6{^-1) + [l-p) 5{^) . (75) 

The reduced flux 4> takes the value of 1 inside the sphere Vs and is equal to outside. 
The probability p that a clump lies inside the domain Ds " from which it may yield a 
signal at the Earth - is just the ratio Ms/JMh of the DM mass Ms confined in that 
sphere with respect to the DM mass Mh contained in the entire domain Vh ■ In the limit 
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where Ad oc rs is small, the DM distribution is homogeneous inside the sphere Vg - with 
constant density Ps{G) - and the probability p may be expressed as the ratio 

^ Mh Mh ^ ' 

For an injected energy Es = 100 GeV and a positron energy at the Earth E = 50 GeV, 
we find a probability p ~ 2 x 10^^ when the statistical domain Vh is chosen to be the 
above-mentioned NFW halo extending up to 20 kpc from the center of the Milky Way. 

Because p is vanishingly small and the number of clumps Nh inside the domain Vh 
exceedingly large, the limit of Poisson statistics is reached. The probability to find n 
clumps inside the sphere Vs is therefore given by 

= ^ exp (-iVs) , (77) 



n 



where Ns = p Nh is the average number of clumps that contribute to the signal 

(n) = Ns = . (78) 

Departures from the statistical law H77|l in the case of a realistic positron propagator will 
be discussed in section Ffi . 1 1 when the number Ng of the clumps involved in the flux at the 
Earth is large whereas the opposite regime will be addressed in section |^J The Poisson 
distribution (|77|l is associated to the variance 

al = {n') = Ns . (79) 

In the hard-sphere approximation, the random part 0^ of the positron flux at the Earth 
- the contribution which the entire constellation of substructures generates ~ is propor- 
tional to the number n of clumps lying inside the sphere Vs- We therefore anticipate 
that the relative variance (Tr/{(t>r) should be equal to the relative variance o'„/(n) of the 
Poisson law H77(l . As a matter of fact, in the limit where Ad is small with respect to 
L - and where the hard-sphere approximation becomes valid - the integrals J7i and Xi 
simplify. If the mass density of reference po is set equal to its solar neighbourhood value 
Ps{&), the ratio Ji/ll boils down to l/Vg so that the exact relation H64() simplifles into 

{^r? fpol! Vsfps(Q) Ns • ^ ' 

We have therefore shown that in the hard-sphere regime, the variance cr^ of the random 
flux (fir is indeed given by the variance ct„ that characterizes the Poisson statistics (|77ll . In 
Fig. 13 the relative variance <Jr/{4'r) - solid curves - and its hard-sphere approximation 
1 / ^/Ns ~ long-dashed lines - are presented together for comparison. Above a positron 
energy at the Earth of 40 GeV, the correct calculation and its hard-sphere limit differ 
by less than ~ 5 x lO^-^. The agreement is remarkable. The diffusion length does not 
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exceed ^1.5 kpc in that case and the hard-sphere approximation succcssfuUy describes 
the statistical properties of the random positron flux (j)^. The relative variance (Jb/B^q 
of the boost factor is also well reproduced by the hard-sphere approximation \/ \/Ns and 
both the short-dashed and long-dashed curves are hardly distinguishable from each other 
at high positron energy E. 

6. A Monte-Carlo approach of the cosmic ray flux fluctuations 

6.1. The large Ns limit and the central limit theorem 

When the average number Ns of clumps that are involved in the signal is large, the 
Poisson statistics (I77II becomes the gaussian distribution 



where S = n — Ns denotes the departure of the number n of substructures inside the 
sphere Ds from its average value Ns. The associated variance is cr„ — \/Ns. We therefore 
anticipate that the flux </>,. will also be randomly distributed according to a gaussian law 
with mean value (0^) and variance 

In order to determine the distribution of probability that drives the random flux 

(f>r - generated by the entire constellation of the clumps lying inside the reservoir Vh - we 
should compute the product of convolution of the Nh distributions of probability V{(p) 
associated each to the flux ip oi a single substructure - or alternatively to its reduced 
flux $ as was discussed in section 14.21 Such a task may seem desperate. However, in 
the large Ns regime, the central limit theorem may be fruitfully applied to solve that 
puzzle. This theorem states that the above-mentioned inextricable product of convolution 
boils down into a gaussian distribution with mean value (0^) = Nh{lp) and variance 
cr^ = Nh {('^3^) — (</')^}- We recognize the expressions (|53f) and (|54|l which have been 
established and numerically computed in section EITI Therefore, the probability to obtain 
a flux (j)r at the Earth may be expressed as 



The probability that the total positron flux (f> at the Earth is enhanced by a factor of B 
with respect to a completely smooth DM distribution is readily obtained as 



where the variance (T_b is given by relation (|52|l . Finally the reduced boost r] = B/B^s 
follows also the same gaussian law 
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Fig. 7. One thousand different realizations of the distribution of DM substructures inside 
the galactic halo has been generated by Monte-Carlo simulation. The injected energy is 
Es ~ 100 GeV. A NFW profile has been assumed with typical scale 25 kpc. A mass 
fraction / = 0.2 is in the form of clumps with mass 10^ Mq. Each histogram corresponds 
to a specific positron energy E at the Earth. The number of realizations - each involving 
271,488 clumps - is plotted as a function of the reduced boost rj = B / B^s- Up to an 
overall factor of a thousand - that corresponds to the number of Monte-Carlo realizations 
- each panel features a numerical estimate of the probability density V{ri). 



with an average value of (rj) = 1 and a variance cr,, = (Tb/Bcs not too different from 
(Jr/ {(j)r) as shown in formula H65(l . 

In order to check our theoretical predictions, we have run a Monte-Carlo simulation 
of the distribution of DM substructures in the Milky Way halo. A thousand different 
realizations have been generated at random assuming a NFW DM galactic halo with a 
fraction / = 0.2 in the form of 10^ M0 clumps. In Fig. |7| the number of realizations is 
plotted as a function of the reduced boost -q for 4 values of the positron energy at the 
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Un ~ ub/ Bcff 


Ns 


50 


20.09 


0.04338 
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65 


20.32 


0.06472 


223.5 


80 


20.48 


0.10966 


78.0 


90 


20.57 


0.19680 


24.2 



Table 1. For each value of the positron energy E at the Earth that has been considered 
in the plots of Fig. [7| we have computed the corresponding effective boost Bcs as well 
as the variance (t,j associated to the reduced boost rj = B / Beg. That variance has been 
derived from relation (|65|l and is in excellent agreement with the rms value of the Monte- 
Carlo simulation. The average number Ns of substructures inside the sphere Vs is also 
indicated. 

Earth. These distributions are the Monte-Carlo counterparts of the gaussian law (|84() 
with a mean value of rj actually very close to 1. In each panel, the rms value of the 
histogram is equal - within a few percent ~ to the variance cr^ = ub/Bcs which we 
have derived from expression (|65() and listed in Tab.^for comparison with the results of 
Fig. For completeness, each histogram has been independently fitted by the Gaussian 
distribution 

-|=»p{-<!y^!} . ,85) 

The amplitude mean value ^ and variance cr arc displayed in each panel and the 
corresponding fitted gaussian is featured by the red curve. The width of each bin is 
Ary ~ 0.04 and since we have generated 10'^ Monte-Carlo realizations, we should obtain 
a value of K = 0.04 x 10^ = 40 for the amplitude. This is quite the case. The mean value 
^ of the gaussian is basically equal to 1 whereas its variance a is very close to the Monte- 
Carlo rms value and to cr^ - see Tab. ^ Because the clumps that are involved in the 
positron signal at the Earth are numerous - the average number Ng is still larger than 
~ 20 even at the highest energy E = %Q GeV - the central limit theorem applies and the 
gaussian distribution (|84|l is an excellent description of the statistical fiuctuations of the 
positron fiux. The question arises now to understand how the distribution of probability 
'P{r]) is modified in the limit where Ns becomes smaller than 1. This is illustrated in 
Fig. El 

6.2. The small Ns regime 

When the diffusion range of positrons is small (for energies close to the emission energy), 
the individual probability distribution 7'i($) = 'P{^) (where $ = (p/tpniax) is strongly 
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peaked at ~ 0. As a result, the probability distribution for the total flux $tot 

generated by the N = clumps of domain Vj^ can be approximated by 

7'iv(*tot) = A^^i($tot) , (86) 

for < $tot < 1 and in the regime where ($tot) is vanishingly small. The proof is 
straightforward. The probability Vn is given by 



7'Ar($tot)=/ ■Pl($)7'Ar-i($tot - $)d$ 

Jo 



(87) 



When 7'Ar_i(<I>tot) behaves qualitatively like 7'i(<i>) and is strongly peaked at a value 
close to 0, two regions dominate the contribution to the integral when $tot < 1, namely 
$ close to (where Vi{^) is large) and $ close to $tot (where VN-i{^tot — *&) is large), 
so that 

^iv($tot) ~ 7'jV-l($tot) + 7'l($tot) (88) 

which proves the property (|86|l . This is illustrated in Fig. |51 where the distributions Vi, 
7^2/2 and 7^3/3 are featured as a function of the total positron flux expressed in units 
of the maximal value i^Smax which a single clump can generate at the Earth. The self 
convolution of Pi{^) has been carried out numerically to yield 7^2 and Vz, assuming 
relation ((HJ- The size of the statistical domain Vh has been fixed by setting a low flux 
cut-off of $mm — 0.001. As is clear in Fig.lHl the three distributions are basically identical 
in the range where $ is smaller than 1. This is so because the probability densities are 
peaked at $ = 0. Should we have chosen a smaller domain Vh and hence a larger value 
for the cut-off the distributions would have been less saturated by their low-flux 

behaviour and relation (|86|l would not have applied. In Fig. 10'^ realizations of a 
clumpy DM halo have been simulated with a substructure mass of 10^ M0. On the 
horizontal axis, the histogram features the boost ratio 77 = B / B^ti which is proportional 
to <i>tot- The resemblance with the analytical distributions of Fig. is striking. The red 
curve which is superimposed on the Monte Carlo results of Fig. ^| corresponds to the 
product 

NhVi{v) = ^ I ps{x)d'x , (89) 

with the same values of / and Mc as in the simulation. On a large portion of the range 
extending from up to _B 11 BcS - therefore for a total flux smaller than ipraax ~ 
relation H86() is a quite good approximation. This regime corresponds to the situation 
where a single clu mp happens to contribute si gnificantly more than the others and is the 



framework of the 



Cumberbatch fc Silkl l|2006(l work. Most of the realizations of Fig. ^| 



correspond to small values of the flux ratio <I>tot • The number of clumps effectively implied 
in the signal is on average very small {Ns <C 1). 
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Fig. 8. Distribution of Bmc/BcS for an energy £' = 80 GeV, given the source energy 
£■5 = 100 GcV, with a mass Mc of individual clumps equal to 10^ Mq, 10^ A/©, 3x10'' Mq 
and 10^ M0. When Mc is small, the clumps are numerous enough for the central limit 
theorem to apply. The mass fraction in clumps is set to the value / = 0.2. The resulting 
distribution is a gaussian, as described in the text. On the other hand, when Mc is large, 
the probability that several clumps contribute to the observed signal is small and the 
observed distribution for Bmc/BcS reflects the one clump distribution -Pi($). 



7. A practical example: HEAT excess 



The putative positron excess reported by the HEAT experiment around ^ 10 GeV is sug- 
gestive of an exotic mechanism such as for example the annihilation of wimps potentially 
concealed inside the halo of the Milky Way. 

The present study basically sketches a well defined frame that can be used in order to 
make predictio ns with respect to available DM candid ates and experimental data. HEAT 
measurements ( Barwick et al.ll997 ; Coutu et al.l2Q01 ) have been widely exploited in con- 



nection to annihilating DM, and we can easily verify how our results translate into phe- 
nomenology. For that purpose, we have chosen an illustration based on a DM candidate 
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Fig. 9. Probability distribution 7'„(<f>)/?7 for n ~ 1, 2 and 3, obtained by consecutive 
convolutions of 7'i((f>). 



that was proposed in the so-called warped GUT theoretical scheme (|Agashe fc Servant 



2004|) . Within such an extra-dimensional modelling, the conservation of a discrete sym- 
metry called Za-symmetry (which is related to the stabilization of the proton, like the 
i?-parity in supersymmetry) , allows the survival of the lightest Za-charged particle, which 
has the properties of a right-handed neutrino called LZP hereafter. This particle is ac- 
tually a Dirac fermion, and given that no matter/anti-matter asymmetry is involved in 
that case, the annihilation rate per volume unit of this species is provided by equation^ 
with (5 = 1/4. We thus have considered a fiducial model in which we fixed the LZP mass 
to 50 GeV, and the extr a-dimension mass paramet er to rriKK — 6 TeV. The cross-section 



formulae can be found in 



Agashe fc ServantI 1)2004 ). and are mainly defined by the isospin 



content of any final state. It is noticeable that about 10% of the annihilation product 
is carried out equally by the three charged lepton/anti-lepton pairs, which can provide 
a relevant contribution to a sharp compone nt close to the wimp mass in the injected 
positron spectrum. We then used the Pythia ( Siostrand et al.ll2001 ) Monte-Carlo to infer 
the positron spectrum associated with all decay and fragmentation processes. 

Besides, in order to calculate the expected positron fraction, we made a further as- 
sumption, which asserts that to each positron generated by wimp annihilation and prop- 
agated to the Earth, an electron is associated with the same spectral information. The 
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Fig. 10. Probability distribution of B/Bcs for an assembly of clumps and the correspond- 
ing distribution for one clump, multiplied by N (see text). The curves are similar in their 
common range, as is expected from an analytical calculation of 7^jv(^') (see Fig. O. 



positron fraction is consequently given by the following expression: 



(90) 



(E) + (tye^^bg{E) + 0e- [E] + <^e- ,bg{E) 

where hg indicates non-exotic contributions (secondary for positrons, both primary 
and secondary for electro ns). For those components, we used the estimates by 
Mnskalenkofc Strongl|l99S). 

We show on Fig. ^2 the results obtained when considering an NFW profile, and 
/ = 20% of the halo mass within a radius of 20 kpc being clumpy. We will not discuss 
how the obtained spectrum is compatible with previous works. Instead, wc want to stress 
here the differences between the naive account for a global and wrong boost factor set 
by the product f x Be, and the correct treatment of the problem that we have presented 
in this paper. To this aim, we used two particular Montc-Carlo simulations with clump 
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masses of IO^Mq and individual boost factor of — 200. This illustrates the situation 
where a small number of clumps contributes, at energies close to E = rnj^zp- In the left 
panel of Fig. ^] the closest clump is found to lie at a distance of 1 kpc to the Earth, 
that corresponds basically to a regime in which B/Bcs < 1 (c/. Fig. I10() . In the right 
panel, the closest clump has a distance to the Earth of ^ 0.25 kpc, which is a much less 
probable configuration, with B / B^s > 5. The corresponding probability is less than 1%, 
as shown by Fig.^| In both panels, the solid blue curves feature the correct treatment of 
the boost factor, while the green lines correspond to the naive shift by a factor of / x B^ 
40 in this case. Notice the discrepancy at low energy in both panels, that is consequent 
to the energy dependence of the correct boost factor Scff. In the right panel, the Monte 
Carlo result is five times larger than the naively boosted flux close to M^zp — 5QGeV, 
as a result of the variance affecting this small Ns configuration. 

This indicates how carefully predictions should be made when computing the flux 
enhancement due to dumpiness. We stress that the variance of that boost should also 
be provided along with the mean values. The spectral distorsions could be sizeable when 
compared to the experimental error bars of the HEAT results. It is of paramount im- 
portance to properly take them into account w hen studying the discover y potential of 



20061) or PAMELA 



the next generation experimen ts, such as AMS (|AMS Collaboration! 
( PAMELA collaborationl200fi K The case of the LZP has been chosen as typical. For par- 



ticles annihilating mostly into charged lepton pairs (respectively quark pairs), - like the 
lightest Kaluza-Klein candidate S'-^-' of universal extra-dimension theories (ncutralinos 
in mSUGRA) - the effect would be stronger (a bit weaker). 

Recently, the pre sence of local dark matter substructure has been invoked with the 



specific assumption (jCumberbatch fc Silk 



2006|) that a single neutralino clump would 



generate alone the observed distortion in the positron spectrum should it be very close 
to the Earth - at a distance of ~ 0.1 pc. The contribution of the other and more remote 
protohalos was assessed to be negligible. 

We will not discuss here whether the distortion which that single substructure gen- 
erates really matches or not the data nor will we be interested in the specific nature of 
the wimp at stake. The real question which we would like to address is to determine the 
probability for a nearby clump alone to shine more strongly than the rest of the other 
protohalos. This situation co uld indeed arise as suggeste d by Fig. ^2 but its probabil- 



ity is vanishingly small in the 



Cumberbatch & Silkl (|200fil ) configuration. The authors of 



that a priori appealing proposal have assumed that half of the Milky Way dark matter 
halo was made of ^ 10^^ E arth mass clumps as suggested by recent numerical simula- 



tions (jPiemand et al 



200,^ . That constellation of neutralino substructures is randomly 
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Fig. 11. Expected positron fraction in the frame of LZP DM particles, with rriLzp = 



50 GeV, superimposed on the HEAT data l|Coutu et al 



200lt) . as computed from two 



random realizations of a clumpy halo with clump masses of IO^Mq. The smooth DM 
distribution follows an NEW profile with a scale radius of 25 kpc, and the mass fraction 
in clumps within 20 kpc has been set to / = 0.2. An intrinsic clump boost of Be ~ 200 
is considered. Left: the closest clump has a distance to the Earth of ^ 1 kpc. Right: the 
closest clump has a distance to the Earth of ~ 0.25 kpc, which is less probable than the 
left-panel case. On both panels, the red and blue curves give respectively the primary 
and the total (primary + background) contributions from DM annihilation. A smooth 
halo would not give rise to any excess with respect to the background in that case. Eor 
comparison, the green curve illustrates the wrong use of the boost, which is simply a 
shift of the spectrum by a factor of / x Be, 40 in that case. 

distributed and contributes on average a positron flux {(pr) at th e Earth. The nearby 



protohalo yields in addition a signal (f>r — {4>r) that, according to ICumberbatch fc Silk 
1 2006(1 ■ overcomes the contribution (0^) from the other clumps. The explanation of the 
HEAT excess in terms of that providential protohalo relies therefore on the assumption 
that the total positron flux (f>r is larger than twice the average value {4>r)- The probability 
for such a conflguration may be expressed as 



PWr>2{(f>r)}= I V( 



+ 00 



(91) 



The distribution of probability is given by the Maxwell law (|82|l insofar as the central 
limit theorem can be applied in that case. The previous relation translates into 



\nP{(l)r > 2{(j)r)} = -a - In 2 - ^ In (tto) 



(92) 



where the parameter a stands for the ratio 
which we need to derive is the variance (t„. 



^2(7^. The relevant statistical quantity 
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Wc may readily apply the tools which wc have constructed up to the condition that 
we are now dealing with a neutralino that produces a continuous positron spectrum 
instead of a line. The positron propagator Ge+ {xq,E x, Es) needs to be replaced by 
the convolution 



G' {x) = Gl+ {xq,E ^x)= Ge+ {xq,E ^ X, Es) 

JE 



dEs , (93) 

Es 



dEe+ 

where dNf,+ /dE^+ denotes the positron spectrum at the source. Without loss of gener- 
ality, we have focused our discussion on a neutralino with mass m-^ = 100 GeV and a 
positron energy at the Earth of = 10 GeV. Neutralinos has been assumed to annihilate 
into bb pairs. According to the hard-sphere approximation of section|Sl the positrons that 
are produced at the energy Es originate from the volume Vs ~ ('\/27r Ad)^ surrounding 
the Earth. In the case which we consider here, the positrons that are detected at the 
energy E at the Earth have been produced at an energy Es that spans the entire range 
from E up to the mass m^. The volume from which the signal originates on average is 
the convolution 



VI, ^ / Vs{E,E, 
Je 



. dN,+ 

s. 



dEs ■ (94) 

Es 



dE^+ 

The number of protohalos that contribute to the positron signal at 10 GeV is inferred to 
be 

and the relative variance (Jr/{4>r) may be crudely approximated by \/ \/Ns- With a 
fraction of clumps of / = 0.5 and a substructure mass of Mc = 10~^ M©, we find a 
number of protohalos of Ns = 2 x lO^'^ and a relative variance of cr ,./((/),.) ^ 2.22 x 10^*". 
Because of that very large number of clumps, the use of the central limit theorem is 
plainly justified. The relative variance of the positron signal is vanishingly small and we 
therefore anticipate that the probability P {<f>r > 2 ((pr)} is completely negligible. 

The correct calculation of the variance makes use of relation (|64(l where the inte- 
grals Xi and J7i are now computed with the convoluted positron propagator G' (a;). We 
have derived a value of <jr/{4>r) ~ 4.13 x 10"^ in good agreement with the hard-sphere 
approximation. We therefore reach the conclusion that 

logioP{(j)r>2{(t)r)} = -1.27 X 10^^ . (96) 

With such an exceedingly small value of the probability, the configuration in which a 
single clump overcomes the signal fro m the other 2 x 10^"^ subst ructures is completely 
unlikely and the hypothesis pursued in 



Cumberbatch fc Silk 



should be abandoned. 



This example illustrates how the tools which we have presented in this article may be 
fruitfully applied in order to derive quantitative results and not just mere qualitative 
arguments. 
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8. Discussion and conclusion 



Summary : The enhancement of indirect signals coming from dark matter annihilation in 
a clumpy halo is usually described by a "boost factor" . We have shown that this quantity 
should be considered as a random variable and we have investigated its statistical prop- 
erties in the following situation: (i) clumps arc distributed like the smooth component; 
(ii) all the clumps arc identical. Wc showed that the boost factor may strongly depend 
upon the specific realization of the clumpy halo we are living in. 



Spatial distribution of clumps: A more realistic model would include the mass profile of 
these clumps, their own density shape and geometry, as well as a number distribution 
inspired by hierarchical structure formation studies (w hich is often found to b e close to 



the smooth distribution, at least at large galactic radii ijBerezinskv et alll200,'j) V Taking 
another number distribution would essentially modify the shape of the effective boost 
factor (the integral Xi of ea. ll9() . whereas individual clump properties would affect mostly 
our estimates of its variance at short scales. For example, the clump number distribution 
is very like ly to be cut off inside the galactic bulge because of strong tidal interactions 
with stars 1 Berezinskv et allboOG ) . However, we do not expect that the results presented 



here would be strongly affected by these effects. For instance, we have shown that the 
steep energy dependence of the mean boost factor, in the case of a positron line, was 
mostly due to contributions from our very local environment, due to the short scale of 
positron propagation. A cut-off in the clump distribution for galactic radii less than ~ 3 
kpc (corresponding to a minimum of 5 kpc from the Earth) would then significantly 
diminish its low energy contribution, and would thus increase the relative variation of 
the boost factor with energy. 



Cumberbatch &Silk: It has been proposed recently (jCumberbatch fc Silkll2006() that the 
positron excess observed by HEAT could be due to the presence of a single clump located 
near the Earth. However, the situation in which the signal due to one clump dominates 
over the background due to all the others is very unlikely. The probability that a clump 
lies in the close proximity of Earth is sizeable only if the density of clumps is high, which 
in turn implies that many of them contribute to the measured flux. We showed that the 
quantitative study of this situation leads to unreasonably small probabilities. 
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